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Abstract. 

An analysis of minijet production from the hard pomeron with a running coupling constant is 
performed. Two supercritical pomerons found in the numerical study are taken into account. The 
calculated inclusive jet production rate is finite at small k± and behaves like 1 /fcj 4 at high k± modulo 
factors coming from logarithmic terms. The average k± is found to be very large (~ 10 — 13 GeV/c) 
and practically independent of energy. This is interpreted as an indication that at present energies we 
are still far from the asymptotics and that, apart from supercritical pomerons, other states contribute 
significantly. 
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1 Introduction. 


In this paper we continue the analysis of a model where a running coupling constant is included into 
the dynamics of reggeized gluons. This model [1,2] is based on the bootstrap relation j|] between the 
reggeized gluon trajectory u>(g) and the interaction kernel for a gluon pair K q (qi , q[ ), which guarantees 
that the production amplitudes satisfy unitarity in the one-reggeized-gluon-exchange approximation 
[Q. To satisfy it both u> and K are written as functionals of a single function 77 (g) which is proportional 
to q 2 in the BFKL fixed coupling case ||. For the running coupling case the asymptotic form 77 (g) « 
g 2 /2a s (g 2 ) at g —> 00 can be derived by comparing the evolution equation with the GLAP equation 
in this limit [ 1 , 2 ]. 

We recall the basic equation for the pomeron with a running coupling: 

(-w(gi) - w(g 2 ))</>(gi) + J ^^K^ ac \ qi , q[)<j>{q[) = E(q)<j>( qi ) ( 1 ) 

where 


w(g) 

K q (qi,q[) 


-a s N c r/(q) 


—2TjT 2 a s 


f d 2 qi 1 
/ (27t) 2 7?(gi)?7(g 2 ) ’ 
/ g(gi) V( Q2) , _ 

vWi) viti) v(qi 


g = gi + g 2 

2 _ g(g) \ 

-gi) »7(gi)»?(g2)/ 


Ti( 2 ) are the colour operators for the two interacting gluons and N c is the number of colours; in the 
vacuum channel XjX^ = —N c . Function <f> is the amputated wave function. We will be interested 
in its computation for the forward scattering (q = 0 ), so dividing by rj once and twice we get the 
semiamputated and the full wave functions, ip and $ respectively. The bootstrap condition is 


/ ^ K i 9luon Hq,Qi,Q[) = u(q)-u(Qi)-u(qi) ( 2 ) 

where in the gluon channel TiT 2 = — N c /2. 

In the following, as in [1,2], we use a parametrization 77 (g) = ( b/2n)(q 2 + TO 2 )ln((g 2 + to 2 )/A 2 ) 
with b = (1/4)(11 — (2/3 )Np) and m > A and we restrict to the physical case N c = 3. We choose 
A = 0.2 GeV for four acting flavours. 

In our work Q we studied the vacuum channel equation (jlj) numerically. Its spectrum relates to 
the pomeron trajectory as a(q) = 1 — E{q) ss 1 + A — a’q 2 . It turned out that two supercritical 
pomeron states appear as normalizable solutions of (|lj). Choosing to = toi = 0.82GeV/ we found for 
them 

A 0 = 0.384, ot’ 0 = 0.250 (GeF/c)" 2 ; Ai = 0.191, a/= 0.124 (Gel^/c )" 2 (3) 

Their wave functions were also computed for the forward case q = 0. 

As an application, total cross sections for the 7 * 7 * scattering, and more qualitatively, for the 
7*73 and pp scattering, were studied in joj] to see the linearization effects. It turned out that the 
contribution from multiple pomeron exchanges became significant only for superhigh energies of the 
order of lOOTelA This seems to support the idea that we are still very far from the asymptopia at 
present energies, so that other states from the two-gluon spectrum give a large contribution. 

In this paper we apply our model with the two found supercritical pomerons to jet production. 
This process has been extensively studied in the framework of the standard BFKL approach with 
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a fixed coupling [7,8]. As well-known, this analysis has lead to some far-reaching conclusions as to 
the importance of mini-jet production at high energies and the logarithmic rise of the multiplicity. 
However in the fixed coupling approach, various ad hoc modifications of the basic BFKL model had 
to be introduced to cut off the spectrum at low k± and also to correctly reproduce the high k± tail of 
the spectrum. Both these problems are naturally resolved by the introduction of a running coupling 
in our model, in which no new parameters appear in contrast to the fixed coupling approach. 

In section 2 we derive the basic equations and study the asymptotic behaviour of the jet production 
cross section. In section 3 we present our numerical results. Discussion and conclusions follow in 
section 4. In the Appendix we study the asymptotic form of the pomeron wave functions. 


2 General formalism. 

In this section we consider the scattering of two highly virtual photons (or heavy ”onia”) to make the 
derivation more rigorous. To obtain the formula for the inclusive jet production let us recall the total 
cross section for one pomeron exchange Q : 

-G(u, 0, qi , q[) I] (! - e iqir )(l - eW) (4) 

i=1,2 

Here q(p) is the momentum of the projectile (target); v = qp = (l/2)s; G(v, q, qi, q[) is the non- 
amputated pomeron Green function in the momentum space, q being the total momentum of the 
two gluons and qi ( q [) being the initial (final) momentum of the first gluon; p q ( p ) is the dipole colour 
density of the projectile (target). For the virtual photons we take the colour dipole densities from 


d 2 rd 2 r'p q (r)p p (r') 


d 2 q\ d 2 ( 


In the high energy limit the dominant term of the Green function comes from the two mentioned 
supercritical pomeron states: 


Gp(v,q,qi,q[) = ^2 v aiM ?$i{qi, &)$*(?£, q 2 ) (5) 

2 = 0,1 

where a* and <1^; are the trajectories and wave functions of the leading (0) and subleading (1) pomerons. 

Minijets appear as intermediate gluon states in the Green function (||). They possess arbitrary k± 
subject to condition k\ « s. To calculate the inclusive cross-section for their production we have 
only to split the Green function in (Q) as indicated in Fig. 1 and drop the integration over k ±, i.e. to 
substitute 


G(v, 0,gi,^)=^ J ^-^^p-G(iy 1: 0,q 1 ,k 1 )V(k 1 ,k' 1 )G{u 2 ,0,k' 1 ,q' 1 )d < - 2 \k 1 +k ± -k[) (6) 

where, for the running coupling case, 

W/7 7/1 c J7(ftl)f7(&l) o /nl 7-71 

v(fc 1 ,«i) = 6 ^_ ~3 r?(0) (7) 

Remembering that qi + q 2 = q[ + q’ 2 = q = 0 we obtain 


I(y,k. l) = 


d 3 


dyd 2 k± 


a = [d 2 rd 2 r' Pq (r) Pp (r') [ e—)(1 - e^) f 


(2 tt) 4 


(2tt)4 


+k ± - k'JGiyu 0, <?!, k\)G(y 2 , 0, k[, q[) - 3r?(0)] 


( 8 ) 
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The appropriate kinematical variables are the rapidity and transverse momentum of the observed 
gluon (jet), y = 5 Inand k±, respectively ; we also have Sj = 2^ and Sip) = k_\_y/se^^ v . It is 
convenient to use a mixed (momentum-coordinate) representation for the Green functions: 

G[y, 0 , r,k) = J j^e^ r G(v, 0 , q u k) (9) 

Defining A$i(r) = <Iq(r) — $^(0) we can write 


We also introduce 


G(u, 0, r, k) — G(v, 0,0, k) = ^ iA A$i(r)$*(fc) 

i= 0,1 


Rf p) = (A®i) q{p) = J d 2 rp q{p) (r)A$i(r) 


( 10 ) 


( 11 ) 


In this notation and also using both the semiamputated and full wave functions we find 

/<»,*x>- £ / |y[ 6 ^ (tl) X) +tx) ~ 3 ’ ,(0) i >, i tl)$,(tl + t - L) ] (12) 


The found supercritical pomerons are isotropic in the transverse space. As a result the inclusive 
cross-section I(y, k±) also turns out to be isotropic. Therefore we can integrate over azymuthal angles 
in (12). Defining for the full and semi-amputated wave functions the integrated quantities 


Hki,k ± ) 



da<&(k\ + kj_) 


i>(k ll k±) = 


darj{ki + k±)$(ki + fcj_) 


(13) 


/ 0 


where a is the angle between k\ and k±, we finally obtain the inclusive cross section 
\ " r -y(A i -A i )( k ± S }h(&i+&i)p'lpP ( 


4 ^ 4 i 0 

i,j=0 


dydk 


dkf 

(27r) 3 . 




V{k±) 


(14) 

The results of numerical calculations of the cross-section (14) and also its generalization to more 
interesting cases of hadronic targets or/and projectiles will be discussed in the next section. In the 
rest of this section we shall study the asymptotical behaviour of the found inclusive jet production 
cross-section at very small and very large transverse momenta and also its y -dependence. 

As to the latter, all y-dependence in ( pT| ) comes from the factor exp(—j/(Aj — A,-)) which has its 
origin in the existence of two different supercritical pomerons. Evidently in the limit s —> 00 this 
dependence dies out, since the relative contribution of the subdominant pomeron becomes neglegible. 
The model thus predicts an asymptotically flat y plateau at very high energies. 

At small k± the cross-section © evidently goes down as k Al , since all other factors are finite in 
this limit. However one should remember that (^j) gives the dominant contribution only while (/;:(_.s) 
continues to be large. At too small k±, when the above quantity becomes finite, all other states from 
the spectrum of the two-gluon equation ([]]), hitherto neglected, begin to give comparable or even 
dominant contribution, so that the found k Al behaviour ceases to be valid. 

To find the asymptotic behaviour of the inclusive cross section for k± —> 00 we need to know the 
behaviour of the pomeron wave functions in the momentum space at q —> 00 and in the ordinary space 
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at r —> 0. In the Appendix we show that 

~ -4 (In V 2 ) 13 

q —»oo 

V’M ~ (ln-)^ 1 ( 15 ) 

r—s-0 x r 

where j3 = — 1 — so that in the forward scattering case j3 depends just on the intercept of the 
corresponding pomeron state. Let us study now the behaviour of (© for k± —> oo. Using $(g) = 
~ ^(l n 9 2 ) /3 1 we find I d 2 q^{q) < oo, so that for the second term in the integrand of ( |l2| ) 

we get 

[ + —■* <S>j(k ± ) [ cPhQiih) ~ -^(Inkl)^ 1 (16) 

J k? ± ^oo J 

To analyze the first term we use the identity 


(2tt) : 


/•OO 

+ k_ l) = 27r / rdrJo(k_ L r)'ip i (r)'i/; j (r) 

Jo 


and the relation 


/*oo ^ n oo 

/ rdrJo(k±r)f(r) = ~— / rdrJi(k±r)f'(r) 
Jo Jo 


which is valid for naturally behaved f, such that 
1 \/3<+/3j+ 2 


rf(r)Ji{k±r) 


= 0 . 


- r—0 


Putting f(r) ~ (in i) * J we obtain the leading behaviour 


cPki-ipifa)ipj{ki +kj_) ~ T( lnfc l) 


/3i+/3j+l 


(17) 

(18) 


(19) 


So at k_ l -toowe find for the inclusive cross-section (12) 

i(y,k±) 


E 

i,j= 0,1 


4 J jfei 


(infel) 




( 20 ) 


This asymptotics corresponds to the standard quark-counting rules behaviour (~ \/k\), modified by a 
power factor due the pomeron energetic dependence and a logarithmic factor coming from the pomeron 
wave function, i.e. from the running of the coupling. Note that in (|20| ) the k± and s dependencies are 
separated. As a result we find that the average (k±) is finite and independent of s and (k 2 ^) ~ s A , 
since it formally diverges for ( 20 ) and one has to restrict k\ < s. 

For the the multiplicity 

(n> = \ / d^ki dyd2k± (21) 

the standard asymptotical behaviour (n) = a Ins + b is obtained. Indeed integrating ( |l2| ) we get 
/ dyd 2 k _l d yji k± = Bs A ° Ins + Cs A °. Since the total cross-section has the form a = As a ° at large s, 
we get the mentioned asymptotical expression for the multiplicity. Of course this result is valid only 
in the extreme limit s —> oo. At large but finite s the existence of two different pomerons leads to 
some additional non-trivial s-dependence. 

It is instructive to compare our asymptotic results with those obtained in the BFKL fixed coupling 
model. In the latter case the inclusive cross section is of course badly behaved in the k± —> 0 limit 
due to scale invariance. It is also very different in the high k± limit. At very large k± such that 
In k± ~ Vln s one finds 


V dydk 2 ± J bfkl 


a(y) 


(k 2 ±S ) A 

k\ 


g- In 2 /° 2 (y) 

In k\ \/ln s 


( 22 ) 
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where a 2 ~ (Ins — 4y 2 /lns). Thus for large k± the BFKL cross section goes down faster than any 
power. Also one obtains that (ln/cj_) ~ Vln s so that both k± and k\ grow with s. 

To bring these predictions in better correspondence with the physical reality, as mentioned, various 
modifications of this orthodox BFKL approach have been introduced. In particular in [7,8] the fusion of 
gluons via ’’fan” diagrams was assumed to take place at high gluonic densities, which was considered as 
a way to partially restore the s-channel unitarity. Then, under some additional assumptions, a 1 /k\ 
asymptotic behaviour similar to (|2(]) was found. In our model such a behaviour naturally follows 
without introducing additional assumptions nor imposing the unitarity restrictions. 

3 Numerical results. 

Taking the wave function evaluated numerically in j^] we have computed the cross section (jl4j). 

In Fig.2 we present c[ 2 a/dydk± for the process 7 * 7 * (in units c = 1). We have chosen the projectile 
photon to have virtuality Q = 5GeV/c and the target one to have virtuality P = IGeV/c. The center 
of mass energy is t/s = 540GeU. 

Of course, processes involving hadronic targets or/and projectiles are much more interesting from 
the practical point of view. However these require some non-perturbative input for the colour densities 
of the colliding hadrons. A possible way to introduce it is evidently to convert Eq. (1) into an evolution 
equation in 1/x and take initial conditions for it from the existing experimental data. Postponing 
this complicated procedure for future studies, we use here a simpler approach, taking for the hadron 
(proton) a Gaussian colour density with a radius corresponding to the observed electromagnetic one. 
Such an approximation, in all probability, somewhat underestimates the coupling of the hadron to the 
pomeron, since the coupling of the latter to constituent gluons is neglected. Nevertheless, we hope 
that it gives a reasonable estimate for the inclusive cross-section. Using this approach we get the 
inclusive jet production cross-sections for the 7 *p and pp scattering shown in Fig.3 and 4 respectively. 

A common feature of jet production in all processes is that the cross-section reaches a maximum 
at k± « 1 GeV/c from which it monotonously goes down both for smaller and larger k±. 

In Fig. 5 the cross-section da/dy integrated over k± in the interval 0.5 -j- 20 GeV/c is presented 
for the process pp. The limitations in the numerical calculations of the wave functions do not allow us 
to study higher values of k± , so that we cannot numerically reach the region where the asymptotical 
behaviour ( | 2 (j| ) is strictly valid. 

In Fig. 6 we show jet multiplicities as a function of s for the three studied processes 7 * 7 *, 7 *p 
and pp. As one observes, their magnitude and behaviour prove to be quite similar. 

We also tried to estimate the average ( k ±). Unfortunately, altough it exists according to (20), 
its value results very sensitive to the high momentum tail of the pomeron wave function, poorly 
determined from our numerical calculations. To avoid this difficulty we chose to calculate the average 
(lnfcj_/A) instead. This average at fixed y depends on s and y only due to the existence of two 
different pomerons, so that this dependence should die out at large enough s. In fact the found 
average (lnfcj_/A) turns our to be practically independent of y and very weakly dependent on s in the 
whole studied range of s and y , rising from 3.97 at y / s = 20 GeV to 4.19 at y 7 s = 20 TeV. These 
values imply a rather high average rising from 10.6 GeV/c at -^/s = 20 GeV to 13.2 GeV/c at 
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Vs = 20 TeV. 


4 Discussion 

Our calculations show that the introduction of a running coupling constant on the basis of the boot¬ 
strap condition cures all the cleseases of the orthodox BFKL approach for jet production. At high fc_L 
the cross-section becomes well-behaved and more or less in accordance with the expectations based 
on the quark counting rules. At small k± no singularity occurs, although contributions from other 
states is expected to dominate. 

As to s- and y -dependence, our predictions are even simpler than in the BFKL approach, since 
the running of the coupling converts the branch point in the complex angular momentum plane, 
corresponding to the BFKL pomeron, into poles, of which only two are located to the right of unity 
and contribute at high energies. As a result at superhigh energies, when only the dominant pomeron 
survives, the y-dependence completely disappears and the s-dependence reduces to the standard s A ° 
factor. At smaller s some y- and extra s-dependence appears due to the existence of two supercritical 
pomerons. 

With all these refinements, some basic predictions of the BFKL theory are reproduced. The 
cross-section for minijet production rises fast and saturates the total cross-section as s —> oo. Jet 
multiplicities rise logarithmically. 

However ( k ±) turns out to depend on s weakly and its calculated value results pretty high, of the 
order of 10-12 GeV/c. This should be contrasted to the experimentally observed much lower values of 
(k_ l) rising with energy. A natural explanation of this discrepancy follows from the conclusion made in 
[6] from the study of the structure functions and total cross-sections in our model: at present energies 
the contribution from the two supercritical pomerons only covers a part of the observed phenomena 
because we are still rather far from the real asymptotics. The bulk of the contribution comes from 
other states, which produce a much softer spectrum of particles. The observed rise of the ( k ±) is then 
related to the dying out of these subasymptotical states. Our prediction is then that the rise of ( k ±) 
should saturate at the level of 10-12 GeV/c. 

This circumstance has also to be taken into account when discussing the numerical results presented 
in Figs. 2-6. Probably they also illustrate predictions for considerably higher energies than the present 
ones, at which, according to the estimate made in [6], they should account for ~ 20% of the observed 
spectra. 

To obtain predictions better suited for present energies one should evidently take into account 
all states present in the spectrum of the pomeron equation (1) and not only the two supercritical 
pomerons. To realize this program an evolution equation in v following from (1) seems to be an 
appropriate tool. As mentioned, it could also effectively take into account the nonperturbative effects 
related to the pomeron coupling to physical hadrons. A work in this direction is now in progress. 


5 Appendix 



5.1 The momentum space behaviour 

The asymtptotic form of the pomeron eigenvalue equation at high momentum is [1,2] 


In ^W(<?) 


am in m 


= eVKl) 


1 


d 2 q'W) 


( 23 ) 


lnm 2 q 2 \nq 2 n J [(q — q') 2 + m 2 } ln[(g — q') 2 + m 2 ] 

where we have put E = (3 /b)e (we consider the case m = mi for simplicity) Let us study the last 
(integral) term, which we denote as D. Changing the variable q' = \q\n we present it in the form 

1 f d 2 Ki/)(\q\K) 


D = 


7r J [(n — k) 2 + m 2 /q 2 ]\n[q 2 (n — n) 2 + m 2 ] 


1 


d 2 Ki/j(\q\(n + k)) 


(24) 


7 r J (n 2 + m 2 /q 2 ) \n(q 2 k 2 + m 2 ) 
with n 2 = 1. In the last form it is evident that the leading terms at q —> oo come from the integration 
region of small n. Then we split the total k space into two parts: k > Kq and k < Kq where k 0 is a 
small number 

K 0 « 1 (25) 

The contributions from these two parts we denote as D\ and Zb, respectively. Our first task is to 
show that Z ?2 cancels the kinetic term in (^) irrespective of the asymptotics of 
With small enough k 0 and for a ’’good enough” xjj 


i/j(\q\{n + k)) ~ 4>(\q\n) = ip(q) 


(26) 


so that we can take it out of the integral over n in Di- Actually (]26|) is our definition of a ’’good” 
function, so that we shall have to check if this condition is indeed satisfied for the found asymptotical 

V’('Z)- 

With ([26]) D 2 simplifies to 


Zb = 


VK?) 


2 2 
r K 0 q 


d 2 q' 


it Jq (q' 2 + m 2 ) \n(q' 2 + m 2 ) 
As q —* oo we assume that also 

Ko q —■<■ oo 


= ip(q) In 


In (/CqQ 2 + m 2 ) 
In m 2 


(27) 


(28) 


Evidently this condition fixes the manner in which kq goes to zero as q turns large. This should be 
taken into account when veryfying condition (|2^). Then (|2?]) gives 

In Kq q 2 


Zb = i>(q) In 


In? 


( 29 ) 


We have also 


But according to 


In In Kq q 2 = ln(ln q 2 + In k,q) 


\nq 2 » | 1 iiKq 


In In K^q 2 = In In q 2 + 


In Kq 
In q 2 


+ 


so that we have 
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Then our final result is 


D 2 = 4>{q) ( ln r ^2 + 0(1/In q 2 ) 

In 


( 30 ) 


The first term exactly cancels the kinetic energy (first) term in the pomeron equation. The correction 
term in (|3(]) evidently is much smaller than eip{q ), since the factor which multiplies the function ip 
goes to zero in (|o|). So we can safely neglect it. 

As a result, the asymptotic equation becomes 

am 2 In to 2 


q 2 In q 2 


= eip(q) + D\ 


with the term D\ given by 


D, = - 


1 f d\'ll>{q')e{{q - q'f - K 2 q 2 ) 


{q - q') 2 ln(g - q') 2 


(31) 


(32) 


where we have omitted the terms m 2 in the denominator because of 
As in [2] we pass to the function x(q) defined by 

X (q) =ip(q)q 2 In q 2 

Multiplying (|3l| ) by q 2 In q 2 we find an equation for \ 

1 [ cPq'xWmq ~ q'Y - «§ q 2 ) 


’him 2 = ix(q) + \ j 


q 2 In q 2 


q' 2 In q' 2 


{q - q') 2 ln(g - q') 2 


(33) 


(34) 


At this point we make a second assumption. Namely we assume that in the integral term of (34) 
values m « q' « q give the dominant contribution (as typical for logarithmic integrals). Of course, 
this assumption is also to be checked for the final asymptotics. With this assumption, we can forget 
about the 9 function and also put the last factor to unity in the integral term of (|34|) . We obtain 


2 , 2 -,^[ r dq' 2 X (q') 

am. In m = ex{q) + / -- 

Jo q'- In q' 


Differentiating with respect to q 2 


dx 

dq 2 


X 


q 2 In q 2 


or 


with a solution 


dx 


d In In q 2 


= -x 


X{q) = A exp ( - 


In In q- 


= A(\nq 2 ) 


^-l/e 


The initial function ip(q) has the asymptotics 


A 


where 


4>(q) = (In q 2 )^ 




(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


Now we have to check that our assumptions are indeed fulfilled for the found asymptotics. 
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Let us begin with the second assumption that the values m « q' « q give the bulk of the 
contribution to the integral in ([34]) . Evidently, in order that the integral be dominated by large values 
of q ', it should diverge as q —> oo. This leads to the condition 

£<0, /?>-! (41) 


So our asymptotics can only be valid for negative energies, that is, for bound states. 

Now for the second part of this assumption. To prove that values q' « q dominate we shall 
calculate the contribution from the region q' » q and show that it is smaller. The corresponding 
integral is 

/•OO 

I = Aq 2 \nq 2 / dq' 2 {In q' 2 ) 13 - 1 /q' 4 (42) 

Jq 2 

In terms of x = In q’ 2 

pOO 

I = q 2 \nq 2 / dxx 13-1 exp(— x) (43) 

J In q 2 


The integral over x can be developed in an asymptotic series in 1/ In q 2 : 

r°° r°° n n „2\/3—1 roo 

/ dxx 13 " 1 exp(— x) = - / a^ -1 dexp(— x) = -=-1- (/3 — 1) / dxx^~ 2 exp(— x) = ... 

J In q 2 Jlnq 2 9 Jlnq 2 

(44) 

^Frorn this we conclude that the integral I has the asymptotics 


I = A{\nq 2 ) p 


(45) 


to be compared to the contribution from the region q' « q which behaves as {lnq 2 )P +1 . We see that 
we have lost one power of In q 2 , so that the region q' » q indeed can be neglected. 

Now to the assumption (|b|). We have explicitly 

4>{\q\{n + k)) = {A/q 2 )(n + K)~ 2 {\\\q 2 {n + k) 2 ) 13 = 


{A/q 2 ){ 1 — 2 n,K — K 2 + 4(nK) 2 )(lng 2 ) /3 (l + (j3{2 tik + K 2 )/hiq 2 ) 


(46) 


I) is satisfied with the behaviour of k as indicated in (28). Indeed take k ~ q 

■ i ^ & v r\11 rtAT'vnnlmn’ i vi (1/1 III) L HTTA 1 Vi v A ^ y ri ^ 


and it is evident that 

with S < 1. Then qn ~ g 1_d —> oo and all the correcting terms in (|46| ) have the order q~ 

Note that this does not mean that (^6|) is quite obvious. It is not valid for, say, the exponential 
function. In fact we have 


exp(ag 2 (n + k) 2 ) = exp(c« 7 2 ) exp(a< 7 2 ( 2 ?zK + k 2 )) 

and since qn is large the second factor cannot be neglected. 

In conclusion, we have verified that our assumtions are fulfilled and therefore the asymptotics (39) 
is correct for negative energies. 


5.2 The coordinate space behaviour 

We are also interested in the coordinate space behaviour of the semiamputated wave function 

V’O) = J ^J 0 {qr)il>{q) (47) 
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To estimate ip(r) for r —> 0 we make use of the asymptotic momentum behaviour found previously 
l|), so we can write 

f qo qdq 


ip(r 


27r 


J 0 (qr)ip(q) 


'qo 


qdq ln(q 2 ) /3 

—*(,r)—— 


(48) 


The first integral for r —> 0 is finite; on the other hand the second integral results to be not 
bounded. Infact using the integration variable y = qr we split the y integration region in two parts, 
q$r < y < yo <C 1 and y > yo] defining the two contribution I\ and I 2 respectively, we have 

rvo 


and 


h~ r~( i-b' 5 ~ o-b" 1 

Jqor V r' 0 V r' 

J [ -Jo(y)(\ny + \n-y ~ (ln-)^ 
J y T ’ r’ 


So the asymptotic small r behaviour will be divergent, precisely 


'ip(r) ~ (in-) 
r—►() v r' 


14 / 3+1 


(49) 


(50) 


(51) 
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7 Figure captions 

Fig.l. The substitution for the Green function necessary to calculate the inclusive jet production 
cross-section. 

Fig.2. Inclusive jet production cross-sections for the process 7*7* at ^/s = 540GeI4 with Q 2 = 
25 (GeV/c) 2 and P 2 = 1 {GeV/c) 2 . 

Fig.3. Inclusive jet production cross-sections for the process 7 *p at = 540 GeV with Q 2 = 
25 (GeV/c) 2 . 

Fig.4. Inclusive jet production cross-sections for the process pp at \fs = 540GeVb 
Fig.5. Cross-sections dcr/dy for the pp process at yfs = 540GeIC 

Fig.6. Multiplicities (n) as a function of the center of mass energy y/s for the processes 7*7* (the 
solid curve),7*p (the dashed curve) and pp (the dotted curve). 
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